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Observations at Very High Energies (VHE, E> 100 GeV) of the BE Lac object Mrk501 taken 
with the High Energy Stereoscopic System (H.E.S.S.) in four distinct periods between 2004 and 
2014 are presented, with focus on the 2012 and 2014 flaring states. The source is detected with 
high significance above ~2TeV in ~13.1 h livetime. The observations comprise low flux states 
and strong flaring events, which in 2014 show a flux level comparable to the 1997 historical 
maximum. Such high flux states enable spectral variability and flux variability studies down to 
a timescale of four minutes in the 2-20 TeV energy range. During the 2014 flare, the source 
is clearly detected in each of these bins. The intrinsic spectrum is well described by a power 
law of index T = 2.15 + 0.06 and does not show curvature in this energy range. Flux dependent 
spectral analyses show a clear harder-when-brighter behaviour. The high flux levels and the high 
sensitivity of H.E.S.S. allow studies in the unprecedented combination of short timescales and an 
energy coverage that extends significantly above 10 TeV. The high energies allow us to probe the 
effect of EBL absorption at low redshifts, jet physics and LIV. The multiwavelength context of 
these VHE observations is presented as well. 
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1. Introduction 

With the exception of two starburst galaxies and a few sources in the Large Magellanic Cloud, 
active galactic nuclei (AGN) are up to now the sole detected extragalactic TeV emitters. Despite 
several classes of AGN existing, the majority of the detections in this energy range consists of 
blazars, in particular the subclass of high frequency peaked BLLac objects (HBL). MrkSOl is 
one of the nearest HBL, with a redshift z ~ 0.034, and was the second extragalactic object to be 
detected as Very High Energy (VHE, E>100GeV) emitter in 1995 [1]. It is known for its strong 
variability at all energies and is referred to as an extreme HBE since its synchrotron peak shifts to 
up to two orders of magnitude higher energies during flares [2]. Thus, a flux dependent spectral 
hardening has been observed in both X-rays and in the GeV-TeV band (e.g. [2, 3, 4]). Because of 
these characteristics, Mrk 501 has been observed extensively at VHE, often during multiwavelength 
(MWE) campaigns to try to characterize its SED with simultaneous data, which are of primary 
importance for such a rapidly variable source. In 1997, it was found in a high state with fluxes 
>1 Crab Unit (c.u.) for several months (e.g. [5, 6, 7, 8, 9, 3]) and the historically highest VHE flux 
(E(> 250GeV) ~ 8.3-10“'° cm“^ s“^) was recorded on April 16, 1997 [3]. Such a prolonged high 
state has not been registered since then, but important flaring activities have been observed on a 
few occasions, like in 2005 [10], 2012 [11] and 2014 (this work). 

2. Observations and Analysis 

H.E.S.S. The H.E.S.S. observations have been caiTied out in four different epoches between 
2004 and 2014. The dataset comprises 34 runs' taken in wobble mode with the HESS I array in 
2004 (4), 2006 (5) and 2012 (4) and with the full HESS II array in 2014 (21). The 2004 and 2006 
runs were part of MWE campaigns and were taken on the nights of June 15-16, 2004 and July 
18-19, 2006. Due to the low state of the source, only upper limits could be derived in the original 
analysis and the results were published separately [12, 13]. The 2012 and 2014 observations were 
triggered by flaring states in the TeV domain detected by MAGIC and EACT, respectively. In 
2012 two runs in two consecutive nights were taken (June 11 and 12), while in 2014 Mrk 501 was 
observed on several nights between June 19-25 and on July 29-30. 

In order to have an homogeneous dataset, only events detected with the CTl-4 telescopes 
were extracted from the 2014 observations, discarding the information coming from CT5. All but 
four runs pass the standard quality cuts, for a total livetime of 13.1 h (12.1 h when corrected for 
acceptance). The mean zenith angle is 63.7° and the mean offset from the pointing position 0.50°. 
Data reduction has been performed using the Model analysis [14] with Loose cuts, chosen to lower 
the energy threshold as much as possible. The Reflected Region Background method [15] was used 
for the background determination of the spectral analysis. An excess of more than 2600 photons is 
detected with high significance. Thanks to the higher sensitivity of the Model analysis, Mrk 501 is 
detected also in the years 2004 and 2006. Spectral analyses have been carried out on different data 
subsets using the forward-folding technique [16] and three different spectral shapes. 

RXTE and Swift-XRT. X-ray coverage has been ensured by RXTE and Swift-XRT. The 
lightcurve of the former has been retrieved from [17], while the data of the latter collected in WT 

*The term "run" refers to a single observation, with a typical exposure of 28 minutes. 
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mode between May 12 and July 30, 2012 (obslDs: 00030793177-00030793201) and between May 
31 and June 28, 2014 (obslDs: 00035023035-00035023047) have been analysed for this work. The 
HEASoft software package v. 6.16 ^ with CALDB v. 20140120 has been used. All the events were 
cleaned and calibrated using the xrtpipeline task and the data in the energy range 0.3-10keV 
with grades 0-2 were analysed. The flux values for the lightcurve (Fig. 1) were calculated from the 
spectra of single snapshots integrating between 2 and lOkeV. These were derived in the following 
way: the data were grouped using the grappha tool to have a minimum of 30 counts/bin and then 
fit using XSPEC v. 12.8.2 with a single power-law model and Galactic hydrogen absorption fixed 
to fiH - 1.58 • 10^° cm“^ [18]. A clear harder-when-brighter behaviour is detected, as also expected 
from the literature (e.g. [2]). Contemporaneous X-rays and H.E.S.S. observations exist for some 
nights, and are typically offset by ~ 90 minutes. On June 12, 2012, Swift-XRT observations start 
immediatly after the H.E.S.S. ones have finished. 

Abastunami Observatory. Optical data are obtained with the 70 cm telescope of the Abas- 
tunami Observatory (Georgia), equipped with an Apogee 6E camera (between 1997-2006 a SBIG 
ST-6 camera). Observations taken with a R Cousins filter cover the 1997-2014 period and have 
been analyzed with the Daophot 11 reduction software using an aperture diameter of 10 . Some 
nights of contemporaneous data with H.E.S.S. exist, and strictly simultanous observations have 
been taken during the flare peak on June 23, 2014. 

3. Results and discussion 

Lightcurves. In the top panel of Fig. 1, the run-wise lightcurve of the 2004-2014 H.E.S.S. 
observations above 2 TeV is shown. Two clear flaiing events can be seen in 2012 and in 2014, along 
with two quiescent states in 2004 and 2006 and a phase of moderate activity in 2014. Mrk501 is 
clearly variable also in X-rays (middle panel) and in optical (lower panel). The comparison of the 
three lightcurves shows no obvious correlation in the flux variations of the different energy bands. 
The optical emission is in an average/low state during both TeV flares and this is true also for the 
strictly simultaneous observations during the peak of the flare in 2014. This can be explained by the 
harder-when-brighter behaviour detected in X-rays, since the variability at lower energy has a lower 
amplitude. On the other hand, the source is in high optical state in 2004 and 2006, when the TeV 
emission is very low. This suggests that at least two zones or mechanisms are necessary to explain 
the flaring events in Mrk501. In X-rays, the brightest event is recorded in the night of June 20-21, 
2014 (MJD 56829), but does not show a correspondent maximum at VHE, as might be expected. 
Given the rapid variability in both energy bands, this could be explained by the non-simultaneity of 
the two observations. The TeV flux is roughly constant during this night and is ~30% higher than 
the one recorded on June 12, 2012. On this date, almost simultaneous X-ray observations show a 
~55% lower X-ray flux. Assuming an SSC emission model, a quadratic correlation between the 
two energy bands is expected. Thus, the source should be able to increase its y-ray emission of 
almost a factor 4 in ~ 90 minutes. This can easily be acieved thanks to a flux doubling timescale 
shorter than 10 minutes [19]. A similar behaviour of high X-ray and low y-ray emission is detected 
on the night of June 22-23, 2014. In this case, though, the TeV flux roughly doubles in the ~ 2 hours 

^http://heasarc.gsfc.nasa.gov/docs/software/lheasoft 
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Figure 1: Top: H.E.S.S. run-wise lightcurve above 2TeV. The flaring states in 2012 and 2014 are 
clearly visible (MJD ~ 56089-90 and 56831-32). It is also apparent that the "low" state in 2014 still 
represent a state of moderate activity when compared to the flux levels in 2004 and 2006. Middle: 
RXTE and Swift-XRT X-ray lightcurves in the energy range 2-lOkeV showing high variability. 
Bottom: optical lightcurve from the Abastunami Observatory. No particularly high flux values or 
strong variability are detected in correspondence of the 2012 and 2014 TeV flares. 


of observations, suggesting a further increase to match the X-ray flux state. The opposite behaviour 
is registered during the TeV flare on June 23-24, 2014 (MJD 56831). The X-ray emission is ~ 20% 
lower than the maximum reached three nights before, while the y-ray flux is ~ 4 times higher. 

More details on the flux variability in 2014 can be found in [19]. 

TeV spectra. The very different flux states and high significance of the detection enable 
a temporally resolved spectral analysis down to very short timescales. Flux dependent spectral 
variability is described in the literature (e.g., [3, 10]), although it has been determined mostly at 
energies below 2 TeV. Analyses of several datasets have been performed (total dataset, single nights 
and single runs during the flares) and the 2012 and 2014 spectra are shown in Fig. 2 {left). Power 
law (PL) spectral fits provide in general poor results and are acceptable only for low-flux nights 
and for some run-wise analysis. This is expected from the curved spectra described in the literature 
(e.g. [3, 8, 10]) and from the consideration that the intrinsic emitted spectrum must have been 
absorbed through the interaction with the extragalactic background light (EBL). Mrk501 spectrum 
extends up to ~20TeV. Despite being a close source, a significant optical depth for EBL absorption 
T > 1 is expected. Hence, although the curved power law (CPE) and exponential cut-off power law 
(ECPL) spectral shapes fit the data satisfactorily, they are not a good description of the intrinsic 
spectra. Eor this reason, in order to determine the source intrinsic spectra in the various flux states, 
fits including EBL absorption have been performed. The model of [20] has been used. The spectra 
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Figure 2: Left: H.E.S.S. observed spectra for the 2014 and 2012 flare and low states. The 1997 
CAT spectra [3] are also shown for comparison. Right: Intrinsic photon indices from EBL absorbed 
PE fits as function of flux. Values are given for single nights (full circles) and for single runs (open 
circles) during the flares. By high fluxes, the spectra are clearly hard. 


are best described by an EBE-absorbed PE. The photon indices are 2.3 ±0.1 and 2.15 ± 0.06 for the 
2012 and 2014 flaring states, and 2.7 ±0.1 for the 2014 low state. The respective normalizations at 
3.5 TeV are 6.7 ±0.3, 15.7 ±0.5 and 2.8 ±0.1 cm“^s“^TeV“^ Eits involving more complex spectral 
models are generally not significantly better anymore. 

The inspection of the residuals of night and run fits shows, however, a possible trend towards 
a spectral softening at the highest energies, especially when lower fluxes are considered. The 
impossibility to detect an actual intrinsic curvature can be due to the too narrow energy range of 
about one decade covered by the spectra. On the other hand, the fact that this trend is mostly 
visible at lower fluxes points towards a harder-when-brighter behaviour also at these high energies. 
Indeed, the plot of the changes of the intrinsic photon indices as function of flux (Eig. 2, right) 
shows that the spectra are clearly hard during the flares. This can be explained as being in the 
vicinity of the maximum of the IC peak, which is known to migrate towards higher energies during 
flares and periods of high activity (e.g. [3, 10]). This is the first time that such behaviour is reported 
for energies extending significantly beyond 10 TeV, up to 20 TeV and for such short time intervals. 
The hard PE intrinsic spectrum during the 2014 flare also indicates that there is no sign of being in 
the Klein-Nishina regime even at these high energies. 

EBL determination. It is possible that the optical depth t derived from [20] for the 2-20 TeV 
energy range at the redshift of Mrk501 (z ^0.034) is not precise. In order to test the model, 
the procedure described in [21] was applied. The sum of the EogEikelihood ratio test^ (EERT) 
profiles peaks at the normalization factor a = 0.93]]Qj4, with a value of ~83, which translate in a 
significance of 9.1 cr. The null hypothesis of a non-existent EBE is hence excluded at more than 9 cr 
level. The derived a is compatible with the model of [20]. Taking into account the contributions 
of the two flares only, a = 0.89];q} 4, which is not compatible with the one calculated in [21] (a = 

^LLRT = -2 • log[£.(Ho)/ where Hq is the null hypothesis, i.e. a = 0, and Ni the hypothesis to be tested, 
a >0. The significance is then calculated as sIlLRT. 
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1.27^q}j). This can be explained by the fact that, altough several sources have been used in [21], 
the major contribution is given by PKS 2155-304 alone. Differences could therefore be explained 
as dependencies on the redshift {zMrksoi ~ 0.034 vs. Zpks 2\55-304- = 0.116) and on the significantly 
different energy range covered (~ 2 - 20 TeV vs. ~ 0.15 - 7 TeV). A further comparison with the low- 
redshift subset composed of Mrk421 and PKS 2005-489 only (<z> = 0.051) is not very meaningful 
because of the large errors on the normalization {aiow-z - [21]) and because the various 

subsets yield results with LLRT values between -2 and 2, hence not significant. On the other hand, 
the present result is perfectly compatible with the one derived in [21] for the high energy dataset 
(^high-energy = which covers the range 0.95 -14TeV. 

LIV studies and QG scale limits. The 2014 flare dafa were exploifed fo search for Lorenfz 
Invariance Violation (LIV) effecls leading fo limifs on fhe Quanfum Gravify (QG) scale. In spife of 
ifs low redshiff (z ~ 0.034), Mrk501 is a very promising source for LIV sfudies due fo fhe hard en¬ 
ergy specfrum. The sfudy of energy dependenf fime-delays wifhin a deferminisfic scenario and wifh 
a Likelihood mefhod was performed following [22]. If consisfs of fhe consfrucfion of a low-energy 
femplafe lighf-curve (2 - 4.5 TeV) injecfed in fhe likelihood til procedure wifh photons of energy 
above 4.5 TeV (~ 500 phofons wifh a negligible background confribufion below 0.5%) and sfudied 
wifh dedicafed simulafions. Two cases were considered: a linear and a quadrafic dependence of fhe 
speed of lighf on fhe phofon energy. As no significanl lime delay was found (more lhan 1 cr from 
zero), limifs on fhe QG scale were computed. The one-sided 95% CL QG limifs including sysfem- 
afic uncerfainfies for fhe sub-luminal (supra-luminal) case are 8.5x10^^ GeV (6.4x10*^ GeV) for 
fhe linear term and 1.15x10'^ GeV (l.OxlO'^ GeV) for fhe quadrafic term. The main uncertainties 
come from fhe femplafe lighf-curve binned til, energy calibralion, knowledge of fhe accepfance 
correcfions and of fhe smearing factors in energy. The fofal confribufion of fhe sysfemalic effecls 
fo fhe error calculalion was conservatively esfimated fo be below 100% of fhe sfafislical one. An- 
ofher fype of systemalics may come from fhe energy-dependenf source effecls which may lead fo 
energy-dependenf lime lags, nof considered in Ibis sfudy. As firsf resulfs, fhese are one of besf con- 
slrainls on fhe linear ferm derived from AGN observations (comparable wifh Ihose oblained wifh 
fhe PKS 2155-304 2006 flare [22]) and fhe besf values on fhe quadrafic term obtained until now 
with GRBs [23] and AGNs [22] data. This is mainly due to the exceptional strength of the Mrk501 
2014 flare leading fo a high number of phofons in fhe energy range of 10-20 TeV. 

4. Conclusion 

H.E.S.S. observations of Mrk501 befween 2004 and 2014 have been presented in a MWL 
confexf. Rapid variabilily in fhe 2-20 TeV energy range during flaring slates in 2012 and 2014 
has been reported. No direcf relalion wifh fhe emission al lower energies has been found. While 
fhe differenl behaviour af X-ray energies can be in principle explained by fhe lime offsel befween 
fhe observafions (Ihis would imply fhe capabilily of fhe source fo increase ifs y-ray flux of a facfor 
~four in ~90 minules), fhe optical emission clearly does nof correlale wifh fhe VHE flux. If is 
in low slale in fhe slriclly simullaneous observations during fhe 2014 flare, while if is significanlly 
higher during fhe TeV low slales in 2004 and 2006. This suggesfs fhe presence of af leasl fwo zones 
or mechanisms able fo produce flaring evenfs in Ihis source. The TeV speclra are hard during fhe 
periods of high flux, and sofler olherwise. In particular, fhe hard simple PL shape extending above 
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20 TeV during the 2014 flare peak indicates the absence of Klein-Nishina suppression at these high 
energies. Studies on the EBL intensity exclude the hypothesis of a non-existent EBL at the 9cr 
level. The normalization factors derived with all datasets (a = 0.93^°-j^) and with the flare datasets 
only {a - 0.89;|;q } 4 ) are consistent with the model of [20]. The quadratic limits on the QM scale 
(1.15x10^' GeV) are the best ones derived from GRB and AGN observations. 
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